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Abstract 

The density of trap states in the bandgap of semiconducting organic single crystals has been 
measured quantitatively and with high energy resolution by means of the experimental method of 
temperature-dependent space-charge-limited-current spectroscopy (TDSCLC). This spectroscopy 
has been applied to study bulk rubrene single crystals, which are shown by this technique to be 
of high chemical and structural quality. A density of deep trap states as low as ~ 10 15 cm -3 is 
measured in the purest crystals, and the exponentially varying shallow trap density near the band 
edge could be identified (one decade in the density of states per ~ 25meV). Furthermore, we have 
induced and spectroscopically identified an oxygen-related sharp hole bulk trap state at 0.27 eV 
above the valence band. 
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I. INTRODUCTION 



The performance of organic thin-film transistors (OTFTs) is steadily improving, and 
the charge carrier mobility, as a key figure of merit, has reached values comparable to 
that of hydrogenated amorphous silicon.-^^ As with most semiconductors, the electrical 
performance is determined to a high degree by modifications of the ideal crystal, such as 
intentional doping and other chemical or structural effects which create electrically active 
states in the band gap. For a thorough understanding of the intrinsic capabilities and 
limitations of organic semiconductors, it is now highly desirable to quantitatively study 
such in-gap states, their density of states (DOS) spectrum, their origin, and their stability. 

Field-effect transistors (FETs) on organic single crystals are well suited for determining 
the surface properties, with mobilities near 20 cm 2 /V s routinely achieved (for a review see 
Gershenson et a/.-). Bulk properties, however, cannot be determined with this surface- 
sensitive method and thus one has to use alternative techniques, such as the time-of-flight 
method as shown in the pioneering work by Karl and co-workers,- thermally stimulated 
current, or space-charge-limited-current (SCLC) measurements. For instance, SCLC mea- 
surements with coplanar electrodes were used by Lang et al— to detect a metastable trap 
state in pentacene single crystals, where the SCLC changes over several orders of magnitude, 
indicating trap filling. 

Additional experimental techniques are used to detect defect states in organic semicon- 
ductors: Recently, oxygen-related states at the surface of naturally oxidized rubrene single 
crystals were detected by photoluminescene measurements.- Kelvin-probe force microscopy 
not only offers a direct imaging of the potential across the channel of an OTFT, but addi- 
tionally allows one to extract the DOS spectrum at the semiconductor/insulator interface.— 

In this study we go beyond the basic concept of SCLC and use the spectroscopic character 
of temperature-dependent SCLC spectroscopy (TDSCLC) as described in the theory by 
Schauer et al. 11 to derive the bulk DOS spectrum in ultrapure rubrene single crystals. We 
found that these crystals of high structural and chemical quality have a broad distribution of 
deep states with a low total density of trap states, in addition to a steep band tail. We also 
use the spectroscopic technique to identify the energetic position of oxygen-induced bulk 
trap states in rubrene single crystals. 
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II. METHOD 



A central aspect of TDSCLC is to exploit the spectroscopic character inherent in the 
temperature dependence of the SCLC due to the energy window associated with the Fermi- 
Dirac statistics.— It is assumed that the SCLC is dominated by the charge carriers that are 
thermally excited from a localized trap into delocalized band states. Therefore the valence 
band edge is seen as separating localized and delocalized states, and it is chosen as the 
reference point for the energy scale (E v = 0), with positive energy toward the midgap. The 
basic equations are Ohm's law in the form j = enQri{{x)F(x) and the Poisson equation 
dF/dx = — en s (x)/(eeo). Here F(x) is the electric field strength along the direction x of 
current flow, j is the current density, ee is the product of the dielectric constant and electric 
permittivity (3.5 for rubrene), /^o is the microscopic band mobility, ri{(x) is the density of 
free carriers, n s (x) is the total density of carriers (free and trapped), which is given by 
the convolution of the density of trap states h(E) in the energy gap with the Fermi-Dirac 
function f(E,Ep,T), i.e., n s = J E h(E)f(E, E-p,T)dE. The shape of the current-voltage 
characteristic j(U) reflects the increment of the space charge with respect to the shift of the 
Fermi energy and thus mirrors the energy dependence of the DOS, 

dn s 1 ee (2m-l) | 

with 

_ B(2m - 1) + B 2 (3m - 2) + d[ln(l + B)]/d\nU 

1 + B{m - 1) ' ^' 

Here L is the thickness of the crystal with electrodes on opposite faces, U the applied voltage, 

k-Q the Boltzmann constant, m = d In j/d\n U the logarithmic slope of the j{U) curve, and B 

contains higher-order derivatives oij(U), B = — [dm/dln U]/[m(m — l)(2m — 1)]. The right 

hand side of Eq. ([1]) can be calculated from the current- voltage characteristic measured at 

only one temperature. For a complete reconstruction of the DOS, however, it is necessary 

to relate a given voltage U to the energy of states which are being filled at this value of U. 

A first starting point is to extract an activation energy Ea from the Arrhenius plot of the 

TDSCLC data for a given U, i.e. = —d\a j / d{k-QT)~ 1 . Because h(E) in general is not 

symmetric around Ep,— the energy Ed of the incremental change of space charge is slightly 

shifted from E&, typically by a fraction of k^T— 

E D = E A+ < 3 ;;■">» k B T. (3) 
(2m — l)[m — l)m 



Here n = —d{E^/k-BT)/d\nU is the derivative of the activation energy with respect to the 
applied voltage. To extract the DOS from the shape of the j(U) curves, it is necessary to 
deconvolute Eq. (jTJ with respect to df/dEp, using a high accuracy deconvolution method 
based on spline functions.— 

III. EXPERIMENTAL DETAILS 

The rubrene crystals were grown by physical vapor transport^ under a stream of ultrapure 
Argon 6.0. The starting material (Aldrich purum) was sublimed three times in vacuum. 
Considerable effort was made to avoid contamination with any organic substances, e.g., 
by using glass tubes cleaned in acids. Typical crystals are platelike, where the direction 
perpendicular to the surface corresponds to the long axis of the orthorhombic unit cell.— 
We note that these high-quality rubrene single crystals show in-plane field-effect mobilities 
at the surface of up to 10 cm 2 /V s. 20 

An optimized sample preparation method was found in a slightly adapted "flip-crystal" 
approach,— benefiting from the minimized sample handling. The thin crystals (preferred 
thickness <2 /mi) are placed on glass substrates with 20 nm Au on 5 nm Cr electrode stripes, 
where they stick by electrostatic adhesion. A 20 nm gold top electrode is then evaporated 
onto the crystal, which is slightly cooled during the evaporation (T mas k = —8 °C) to minimize 
thermal damage. The overlap of the electrodes results in a typical measurement cross section 
of A ~ 2 • 10~ 5 cm 2 . Finally, electrical connections to a sample holder were made with silver 
epoxy and 25-/mi-thick gold wires. The thickness of the crystals was measured by atomic 
force microscopy, as optical inspection turned out to be unreliable for the ultrathin crystals. 

The electrical measurements were performed in a closed-cycle cryostat in inert helium 
atmosphere and in darkness, covering a maximal temperature range of 30-350 K. By the use 
of helium exchange gas, a constant and reliable sample temperature, measured at the back 
side of the sample holder, could be adjusted with a resistive heater coil (driven by a Lakeshore 
331 controller). A Keithley 6517A electrometer was used for the electrical measurements. 
By a proper shielding, leakage currents well below lpA at 100 V were achieved. To inject 
holes from the bottom contact, a negative voltage was applied at the top contact; V > 
results in hole injection from the top contact. The superiority of laminated vs evaporated 
contacts was reported previously.— 1 ^ 
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The typical measurement procedure for TDSCLC was as follows. First, initial tests at 
room temperature were performed in order to check the reproducibility of repeated mea- 
surements of the j(U) characteristic. Thereafter the sample was cooled to typically 100 K 
at a rate of 2K/min. Prior to the measurement of the current- voltage characteristic (at 
each chosen temperature), an initial delay of 20min ensured thermal equilibrium. The 
j(U) curves were measured by a stepwise increase of the applied voltage and measuring the 
quasiequilibrium current. A current measurement delay of 10 s turned out to be sufficiently 
long compared to the settling time of the system. Typically 50 points were measured per 
voltage decade. The maximal voltage was limited such that the current density would not 
exceed 0.1 A/cm 2 in order to avoid crystal damage. 

Usually j(U) was measured every 10 K between 100 and 200 K. At lower temperatures 
it is difficult to get accurate j(U) curves, because the current increases extremely rapidly 
with voltage. At higher temperatures, the broadening of the Fermi statistics has an adverse 
influence on the spectroscopic character of TDSCLC. In particular, shallow states can only 
be reliably measured at low temperature. 

The issue of possible long-term charge trapping as result of a j(U) sweep needs careful 
attention during the course of a measurement. In most crystals, subsequent sweeps at the 
same low temperature yield identical curves. This indicates a detrapping of the charge on 
a time scale faster than that of a sweep (~ 10 min). For a few samples, however, it was 
necessary to warm the crystal to room temperature after every sweep in order to restore the 
initial condition. Subsequent j(U) sweeps at the same low temperature then yield the same 
results. If charge trapping or sample deterioration influence the measurement, the current 
at a given fixed voltage will not be thermally activated (as opposed to the Arrhenius plot in 
Fig. [2]). Therefore, the measurement itself is a valuable self-consistency test. 
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FIG. 1: Space-charge-limited-current (SCLC) density vs applied voltage at different temperatures 
for a rubrene single crystal (Ru65-2, L = 0.6 /xm). The temperature step is 10 K. The inset 
shows j(U) at 300 K. The straight line indicates the Ohmic behavior of thermally generated charge 
carriers. 

IV. EXPERIMENTAL RESULTS 
A. DOS in rubrene single crystals 

In Fig. (U the j(U) curves at different temperatures are shown, measured perpendicular 
to the molecule layers in a rubrene single crystal (Ru65-2, L = 0.6 /im). The "Ohmic" 
region at low voltages is indicated by a straight line (see inset of Fig. [T]). The onset of SCLC 
at ~ 0.1 V indicates that holes are effectively injected from the laminated gold contact. At 
lower temperatures the number of thermally excited carriers decreases exponentially and the 
Ohmic region disappears below the sensitivity of the measurement setup (~ 10~ 14 A). 

In Fig. [2] the Arrhenius plots of the current density are shown for the same data set as 
in Fig. dJ For clarity we show only data points for selected voltages. The slope yields the 
activation energy E\(U) for each applied voltage U and we note the high quality of these 
plots. The resulting Ea{U) is given in the inset and is in agreement with recent FET data—. 
The effective Fermi energy is moved from ~ 0.45 eV at the lowest voltage to ~ 0.1 eV at the 
highest injection voltage. The distinct change of slope near 0.27 eV is worth noting, as it 
obviously reflects a marked increase of the trap density. We note that, among other evidence, 
the smooth variation of E\ at low U indicates space-charge-limited-transport rather than 
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FIG. 2: Current density for a fixed applied voltage U vs the inverse temperature for the data 
plotted in Fig. [TJ The straight lines are the Arrhenius fits used to determine the activation energy 
Ea(U). The inset shows the resulting E\(U). 

transport limited by contacts. 25 

The analysis described above is used to extract the density of trap states from these 
TDSCLC measurements. In order to calculate higher-order derivatives of the experimental 
j(U) and E\(U) curves according to Eq. (pQ), a smoothing spline fit was applied to the 
measured data, keeping the fit within 1% of the raw data. The resulting density of trap 
states after the deconvolution is shown in Fig. [31 The edge of the valence band is used as 
the reference level and the positive energy axis points to the center of the band gap. Three 
main features can be discerned. (1) An exponential increase of the DOS toward the band is 
observed for all rubrene single crystals. However, the characteristic energy k-gTt over which 
the DOS is reduced by a factor of e varies from sample to sample. The sample Ru52-3 
has a broad exponential distribution with ksT t = 2 10 me V and the highest density of trap 
states. Broad tail states with k-gT t = 180 meV were recently reported for pentacene single 
crystals 2 ^. The sample Ru65-1 with the lowest deep trap density of order 10 15 cm -1 eV -1 in 
the energy range from 0.45 to 0.1 eV has a characteristic exponential distribution parameter 
of ksTt = 180 meV. (2) Of particular interest in this sample are the shallow trap states 
below ~ 0.1 eV with a steeper slope (fce^t = llmeV), reminiscent of band tail states. Due 
to the large increase of the trap density the quasi Fermi level is pinned close to the band 
edge, and the amount of injected charge cannot fill these shallow tail states. In this energy 
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FIG. 3: (Color online) Density of states above the valence band for four rubrene single crystals. 

range the charge transport is still activated ~ 0.05 eV). (3) The fine structure of the 
DOS indicates small features for all crystals, due to discrete trap levels in the band gap. 
The observed fine structure was confirmed using the energy dependence of the statistical 
shift. 11 

In comparison, data of organic thin films show a very similar qualitative DOS spectrum,— 
exhibiting tail states and exponentially increasing deep states, albeit at much higher con 
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centrations than in the purest rubrene single crystals. Again, the FET geometry in Ref. 
might have emphasized defects near the interface, because FET devices made from the same 
type of high quality rubrene single crystals show approximately three orders of magnitude 
higher interface trap density than the bulk trap density presented here.— 

It is remarkable that the representative selection of rubrene crystals vary in their DOS, 
although grown under basically identical conditions.— This variation may originate from 
individual micro-conditions during and after crystal growth, i.e., different actual growth 
temperature, growth rate, (thermo)mechanical strain, and different atmospheric conditions 
in the device fabrication process. On the other hand, measurements of different cross sec- 
tions on the same crystal result in virtually identical DOS spectra, which is a convincing 
verification of the data evaluation. 
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FIG. 4: (Color online) Density of trap states in rubrene single crystals before and after exposure 
to 1 C , 2- The oxygen-induced defect acts as a hole trap at an energy of 0.27 eV above the valence 
band edge. Sample Ru65-2 was illuminated in an oxygen atmosphere to form the 1 C , 2 directly at 
the rubrene crystal surface. Sample Ru71-4 was exposed to oxygen excited by uv light. 

B. Oxygen-induced bulk trap states 

In order to understand the role of bulk traps in organic single crystals and to demonstrate 
the power of TDSCLC spectroscopy, we investigated the influence of oxygen on the trap 
density of rubrene single crystals. It is well known that the reaction of rubrene with singlet 
oxygen 1 C>2 forms the endoperoxide.— Rubrene itself can act as the sensitizer when the triplet 
state of rubrene is populated by an intersystem crossing from the singlet state excited with 
visible light. The energy of this triplet state at 1.2 eV above the ground state is transferred 
to the molecular oxygen resulting in 1 2 which reacts with the rubrene molecule. 

After measuring the full DOS in the as-grown crystals (open symbols in Fig. Hj), we 
illuminated the sample Ru65-2 with visible light under oxygen atmosphere for four hours. 
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For comparison the sample Ru71-4 was directly exposed for four hours to 1 C>2 by exciting 
molecular oxygen with uv light in the vicinity of the sample; the sample itself was held 
in the dark. The exposure of rubrene to 1 C>2 results in a large peak in the density of 
trap states at 0.27 eV above the valence band (filled symbols in Fig. H]). The amount of 
oxygen-induced trap states in sample Ru65-2 is estimated to be N° x » 2 • 10 17 cm~ 3 , which 
corresponds to ~ 100 ppm. Because the sample was illuminated in an oxygen atmosphere 
without removing from the cryostat we can exclude origins other than oxygen for this hole 
trap. We note that the energetic position of the identified oxygen trap state is in agreement 
with recently published photoluminescence measurements on oxydized rubrene.- In sample 
Ru71-4 the oxygen- induced trap states are within the measurement error at the same energy, 
but the density of the trap states is slightly higher than for the illuminated sample (N° x w 
3 • 10 17 cm -3 ). In addition a shoulder appears in the energetic distribution on the side closer 
to the valence band. This difference may be due to the formation of O3 under uv light which 
might also react with the rubrene single crystal. 

The trap level at E = 0.27 eV already exists in the pristine rubrene single crystals, with 
much lower concentrations. This is due to the device fabrication process; the samples are 
handled in room air under microscope illumination, which results in a similar reaction as 
for the illuminated sample Ru65-2. The oxygen-related trap level was also observed in other 
rubrene samples with various concentrations, probably as a result of a different exposure 
time to air during handling. 

V. CONCLUSIONS 

In conclusion, we have implemented temperature-dependent space-charge-limited-current 
spectroscopy and demonstrate it to be a powerful tool to quantitatively measure the density 
of bulk trap states with high energy resolution. Applied to high-quality rubrene crystals 
this method reveals the existence of states within ~ 0.1 eV of the band edge, reminiscent of 
band tails, and a smooth distribution of states deeper in the gap. Discrete peaks are also 
observed, and through a controlled exposure of the crystals to activated oxygen, a distinct 
and stable trap level at 0.27 eV has been created. Applying this approach to other organic 
semiconductors will be very helpful in the quest to identify intrinsic and extrinsic factors 
that dominate charge transport in organic semiconductors. 
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